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. We analyze families of non-autonomous systems of first-order ordinary differen- 

tial equations admitting a common time- dependent superposition rule, i.e., a time- 
dependent map expressing any solution of each of these systems in terms of a generic 
set of particular solutions of the system and some constants. We next study relations 
of these families, called Lie families, with the theory of Lie and quasi-Lie systems 
and apply our theory to provide common time-dependent superposition rules for 
certain Lie families. 
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1 Introduction. 

The theory of Lie systems [H |2l [3l HI |5l El [7] deals with non-autonomous systems of first- 
order ordinary differential equations such that all their solutions can be written in terms of 
generic sets of particular solutions and some constants, by means of a time-independent 
function. Such functions are called superposition rules and the systems admitting this 
mathematical property are called Lie systems. Lie succeeded in characterizing systems 
admitting a superposition rule. His result, known now as Lie Theorem states that a 
non-autonomous system (time-dependent vector field) Xf is a Lie system if and only if 
there exists a finite-dimensional Lie algebra of vector fields Vq such that Xt G Vq for all t. 

Note that a superposition rule can be found explicitly even for systems whose general 
solution is not known, like in the case of Riccati equations [8] , and its knowledge enables 
us to obtain the general solution out of certain sets of particular solutions in an easier 
way than solving directly the system. 
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In the theory of Lie systems various methods have been developed to obtain super- 
position rules, time-dependent and time-independent constants of the motion, exact so- 
lutions, integrability conditions, and other interesting properties for particular systems 
P [TOl im, [121 HSl HH [15]. Unfortunately, being a Lie system is rather exceptional and, 
in order to apply the methods of the theory of Lie systems to a broader set of non- 
autonomous systems, some generalizations of this theory have been proposed. The gener- 
alized methods are presently used to investigate some partial differential equations [7], a 
class of second-order differential equations (the so-called SODE Lie systems [lO]), certain 
Schrodinger equations [16], etc. 

With the same aim of applying the theory of Lie systems to a broader family of 
systems, it has been recently developed the theory of quasi-Lie schemes and quasi-Lie 
systems |T71 [TH [19]. This theory allows us to investigate some non-Lie systems and it can 
be applied to dealing with certain second- and even higher-order systems of differential 
equations. For example, it enables us to analyze some non-linear oscillators [17], dissi- 
pative Milne-Pinney equations [IB], Emden- Fowler equations [IH], etc. One of the main 
results obtained through quasi-Lie scheme approach is the existence of the so-called time- 
dependent superposition rules, that is, time-dependent superposition functions expressing 
the general solution in terms of a generic family of particular solutions of this system. 

Note however that the concept of time-dependent superposition rule does not make 
much sense for a single non-autonomous system. This is because, as explained in jl7j . 
any single non-autonomous system admits such a superposition rule which, however, can 
be as difficult to finding as the general solution of the system and, therefore, it cannot 
be generally used to analyze properties of the system. This is analogous to the fact 
that each autonomous system is automatically a Lie system, as each single vector field 
spans a one-dimensional Lie algebra. Therefore, it only makes non-trivial sense to speak 
about common time- dependent superposition rules for a bigger family of non-autonomous 
systems. 

In this paper, we give a natural generalization of Lie Theorem characterizing Lie 
systems. This enables us to show that many families of non-autonomous systems of 
first-order ordinary differential equations are Lie families, that is they admit common 
time-dependent superposition rules. Furthermore, we study some Lie families and we 
obtain common time- dependent superposition rules for all of them. 

The organization of the paper goes as follows. Section 2 describes common time- 
dependent superposition rules in terms of certain horizontal foliations. In Section 3 
we generalize Lie Theorem to characterize families of systems admitting common time- 
dependent superposition rules. We posteriorly use this result to analyze the relations 
between quasi-Lie systems. Lie systems and time-dependent superposition rules in Sec- 
tion 4. We finally apply all our results to investigate some Lie families throughout Section 
5. 
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2 Time-dependent superpositions and foliations. 



In this Section we develop the concept of common time-dependent superposition rule for a 
family of non-autonomous systems of first-order ordinary differential equations and relate 
this concept to certain horizontal foliations. For the sake of simplicity, we investigate 
these concepts in local coordinates, but our approach can be slightly modified to handle 
systems on manifolds. 

Consider a family, parametrized by elements a of a set A, of non-autonomous systems 
of first-order ordinary differential equations on of the form 

— = Y^{t,x), i = l,...,n, aeA. (1) 

In applications, A is often a finite subset of M or A = C°°(R). Solutions of these systems 
are integral curves of the family {Ya}aeA of time-dependent vector fields on M" given by 

" d 

Y^{t,x) = Y,Y^{t,x) — , aeA. (2) 

i=l 

Note 1. In order to simplify the terminology, we will use Yq, to designate both: a time- 
dependent vector field of the above family and the non-autonomous system describing its 
integral curves. 

Denote with the autonomization of the time-dependent vector field Y^, that is, the 
vector field on M x M" defined by 

F.(«..) = | + X:y-(,x)^. 

i=l 

Integral curves of ([1]) can be identified with trajectories of the vector field (autonomous 
system) Y^. Let us state the fundamental concept studied throughout the paper. 

Definition 2. We say that the family of non- autonomous systems ([1]) admits a common 
time- dependent superposition rule, if there exists a map $ : M x — )■ M", 

X = . . . , fci, . . . , kn), (3) 

such that the general solution x{t) of any system Yq, of the family ([1]) can be written, at 
least for sufficiently small t, as 

x{t) = <l>(t,X(i)(t),...,X(„)(t);A;i,...,fc„), 

with {x(a) (t) I a = 1, . . . , m} being a generic set of particular solutions of Y^, and fci, . . . , fc„ 
being constants associated with each particular solution. A family of systems ([T]) admit- 
ting a common time-dependent superposition rule is called a Lie family. 
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Note 3. We do not want to formalize precisely what 'generic' means in the above def- 
inition, as it is not crucial for our purposes and depends on the context. One can have 
in mind the following example: for a system of linear homogeneous differential equations 
'generic' means that the particular solutions are linearly independent. 

Given a common time-dependent superposition rule $ : R x ]R"-(™+i) — > of a Lie 
family {Y^jaeA, the map <l>(t, . . . , ■) : M" — > M", X(o) = . . . , /c), 

is regular for a generic point . . . € M x M"™ and, in view of the Implicit 

Function Theorem, it can be inverted to write 

k = ^(t,X(o), . . . 

for a certain map : M x M^^^+i) — ). M", and k = (ki, . . . , kn) being the only point in 
such that 

X{o) = • • • ,a^{m); k). 

Note 4. As a matter of fact, the maps \E' and $ are defined only locally on open subsets 
of M X but, for simplicity, we will write M x M"('"+i) for their domains. 

Consequently, the map determines locally a n-codimensional foliation ^ of the man- 
ifold M X M"("'+^) into the level sets of ^. Moreover, as the fundamental property of the 
map establishes that ^(t, a;(o)(t), . . . , X(^m)it)) is constant for any (m + l)-tuple of par- 
ticular solutions of any system of the family ([1]), the foliation determined by \E' is invariant 
under the permutation of its (m + 1) arguments {x(^a) | a = 0, . . . , m}, and, differentiating 
\E'(t, X(o)(t), . . . with respect to t, we get 

^ + EE^"(^'^WW)^ = 0' J = l,...,n, aeA, (4) 

a=0 i=l {«) 

where \E' = (\E'\ . . . , \E^"). 

Definition 5. Given a time-dependent vector field Y = 'Y^^-iY^{t,x)d / dx"^ on M", we 
define its prolongation to M x as the vector field on M x M"-(™+i) given by 

m n Q 

F(t,X(o),...,X(™)) = ^^y\t.x^a))^r--. (5) 

a=0 i=l ^^('^) 

and its time-prolongation to M x M"(™'+i) as the vector field on M x M"('"+^) of the form 

n m n 



The equalities (j4]) show that the functions {\E'* | z = 1, . . . , n} are first-integrals for the 
vector fields {l^jagA, that is, Y^^^ = for i = 1, . . . , n and a G A. Therefore, the vector 
fields Yq, are tangent to the leaves of ^. 
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The foliation ^ has another important property. If the leaf Jfc is the level set of 
corresponding to a certain k = G M", and given {t,X(i), . . . ,X(rn)) G M x 

M™", there is only one point a;(o) G such that (t, a;(o), a;(i), . . . G ^fc- Thus, the 

projection onto the last m ■ n coordinates and the time 

tt: (t,X(o),...,a;(™)) GMxM"^'"^!) ^ (t, . . . , X(„)) G M x R"™, 

induces a local diffeomorphism from the leaf of ^ into M x R"™. We will say that the 
foliation ^ is horizontal with respect to the projection vr. 

On the other hand, the horizontal foliation defines the common time-dependent su- 
perposition rule without referring to the map \E'. Indeed, if we take a point X(o) and 
m particular solutions, . . . , for a system of the family, then X(o){t) is the 

unique curve in R" such that the points of the curve 

(t , a;(o) (t) , (t) , . . . , (t) ) C R x M"'" 

belong to the same leaf as the point (0, a;(o)(0), X(i)(0), . . . ,a;(m)(0)). Thus, it is only the 
horizontal foliation ^ that really matters when the common time-dependent superposition 
rule is concerned. It is in a sense obvious, as composing \I' with a diffeomorphism on 
R" changes the superposition function (rearranges the level sets) but yields the same 
superposition rule. This proves the following (cf. [7]). 

Proposition 6. Giving a common time- dependent superposition rule ^ for a Lie family 
(Oy is equivalent to giving a foliation which is horizontal with respect to the projection 
TT : R X R('"+i)" — ). R X R"™ and such that the vector fields {Yaja^A are tangent to their 
leaves. 

3 Generalized Lie Theorem. 

It is generally difficult to determine whether a family ([8]) admits a common time-dependent 
superposition rule by means of Proposition O It is therefore interesting to find a charac- 
terization of Lie families by means of a more convenient criterion, e.g. through an easily 
verifiable condition based on the properties of the time-dependent vector fields {Ya}aGA- 
Finding such a criterion is the main result of this section. It is formulated as Generalized 
Lie Theorem. 

We start with three lemmata. The proofs of first two of them are straightforward. 

Lemma 7. Given two time- dependent vector fields X and Y on R", the commutator 
[X, Y] on R X R"('"+i) is the prolongation of a time- dependent vector field Z on R", 
[X,Y] = Z. 

Lemma 8. Given a family of time- dependent vector fields, on R", their 

autonomizations satisfy the relations 

r 

[Xj, Xk]{t, x) = ^ fjki{t)Xi{t, x), j. A; = 1, . . . , r, 
1=1 
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for some time- dependent functions fjki : R — >■ if and only if their time-prolongations 
to ]R X Xi, . . . , Xr, satisfy analogous relations 

r 

[Xj, Xk]{t, x)^Yl fMt)Mt, j,k^l,...,r. 

1=1 

Moreover, X][=i fjki{t) ^ for all j,k ^ 1, . . . ,r. 

Lemma 9. Consider a farnily of time- dependent vector fields, Yi, . Yr, with time pro- 
longations to R X Yi, . . . ,Yr, such that their projections 7r*(lj) are linearly inde- 
pendent at a generic point inRx W^. Then, JJj=i bjYj, with bj e C~(K x W""), is of the 
form Y (resp. Y) for a time- dependent vector field Y on W^, if and only if the functions 
bj depend on the time only, that is, bj — bj{t), and Y^^j=i — (resp., Yl^j=i — 

Proof. We shall only detail the proof of the above claim for bjYj = Y, as the proof 
of the other case is completely analogous. Let us write in coordinates 

5 " d 

1=1 a=0 (") 

Then, 

r r n m r> 

O 



j=l j=l i=l a=0 ^^{a) 



+ ^6j(i,a;(o),...,a;(„)) — , 



which is a prolongation if and only if there are functions : R x — > R, with i — 
1, . . . , n, such that 

r 

X] bj {t, a;(o) , • • • , x^^))A] (t, ) = B\t, a;(„) ), 

\ a = 0,...,m, i = l,...,n. 

X]6j(i,a;(o),...,X(^)) = 0, 

If the functions bi, . . . ,br are time-dependent only and Yl^j=i — 0; above conditions 
hold and Y7j=i bjYj is the prolongation to R x R'*("*+^) of the time-dependent vector field 
Y ^Y^'^^^B'{t,x)d/dx\ 

Conversely, suppose that Xlj=i ^j^j ^ prolongation for a time-dependent vector field 
on R**. In this case, the functions bj{t, X(o), . . . , X(^rn)) solve the following system of linear 
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equations in the unknown variables Ua- 

( ^ 

r 

J^Uj = 0, 
I j=i 

where a = 1, . . . ,m and i = 1, . . . ,n. This is a system of m ■ + 1 equations and, as 
7T^{Yj), with j = 1, . . . ,r, are hnearly independent by assumption, the solutions are 
uniquely determined by the variables {t, . . . , and therefore they do not depend 
on X(o). Since time-prolongations are invariant with respect to the symmetry group Sm+i 
acting on ]R"(™-+^) = (M'^)™"'"^ in the obvious way, the functions X(o), • . . with 
j = 1, . . . ,r, must satisfy such a symmetry. Hence, as they do not depend on X(o), they 
cannot depend on the variables , . . . , } and they are functions depending on the 
time only. □ 

Theorem 10. (Generalized Lie Theorem) The family of systems ^\) admits a com- 
mon time- dependent superposition rule if and only if the vector fields {FajoGA can be 
written in the form 

r 

Y^{t,x) = Y,bajit)Xj{t,x), aeA, (6) 
i=i 

where b^j are functions of the time only, Yl^=i^oij = 1; o.nd, Xi, . . . , X^, are time- 
dependent vector fields such that 

r 

[Xj, Xk]{t, x) = J2 fjki{t)Xi{t, x), j,k = l,...,r. (7) 

1=1 

for some functions fjM : M — > M, with j, k,l = 1, . . . ,r . We call the family of autonomiza- 
tions, Xi, . . . , Xr, a system of generators of tlie Lie family. 

Proof. Suppose first that the family of systems ([T]) admits a common time- dependent su- 
perposition rule and let be the corresponding n-co dimensional horizontal foliation. The 
vector fields {Ya}a<^A are tangent to the leaves of the foliation ^ and span a distribution 
"Do on M X M"('"+i). Such a distribution need not be involutive, see examples in Section 
|5l Nevertheless, we can enlarge the family {Ya}a<^A to the Lie algebra of vector fields 
generated by such a family. This Lie algebra is spanned by {YajaeA and all their possible 
Lie brackets, i.e., 

y;, [Y^.Ypl [Y^,[Yp,Y,]l [Y^,[Yp,[Y,,Ys]]i... a, (3,j, S, . . . E A. (8) 

All the above vector fields are tangent to the leaves of the foliation ^ and therefore there 
are up to m-n+1 linearly independent ones at a generic point of M X M'^("^+i) . Consequently, 
they span an involutive generalized distribution V with leaves of dimension r < m-n + 1. 
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In a neighborhood of a regular point of this fohation, take now a finite basis of vector 
fields from the elements of the family ([8]) spanning the distribution. By construction, at 
least one of them must be of the form Xi for a certain time-dependent vector field Xi 
on and, in view of lemma |9] and the form of the family ([8]), those not being time- 
prolongations are just prolongations. Therefore, if we add Xi to those elements of the 
basis being prolongations, we get a new basis of the distribution V made up by certain 
r < m-n+1 time-prolongations Xi, . . . , Xr- In other words, the distribution V is locally 
spanned, near regular points, by time-prolongations, say Xi, . . . ,Xr. As the generalized 
distribution D is involutive, there exist real functions fjku with j, fc, / = 1, . . . , r, on 
M X R"('^+i) such that 

r 

[Xj,Xk]=^fjkiXu j, A; = l,...,r, 

1=1 

and as the left side of the above equalities are prolongations, we get that, in view of 
lemma [9], all the functions fjki depend on time only and Yll=i fjki = 0. Finally, taking 
into account Lemma [H we have 

r 

[Xj,Xk]{t,x) = ^fjki{t)Xi{t,x), j, /c = l,...,r. 
1=1 

Note that as the vector fields {Ya}a€A are contained in the distribution V, there exist 
some functions baj G C°°{W^^"^~^^'^) such that = Y7j=ibajXj for every a E A. In 
consequence, according again to lemma [HI the functions b^j depend on the time only, i.e., 
baj = baj(t). Therefore, we get that 

r r 

Let us prove the converse. Assume that we can write 

r 

for certain time-dependent vector fields Xi, . . . , on such that 

r 

[Xj, Xk]{t, x) = Y^ fjki{t)Xi{t, x), j, /c = 1, . . . , r. 

1=1 

In view of lemma [H the vector fields Xi , . . . , X,. span an involutive distribution V on 
M X M"(™+i) for any m. Furthermore, the rank of this distribution is not greater than 
r and therefore, for m big enough, this distribution is at least n-co dimensional and it 
gives rise to a foliation which is horizontal with respect to the projection vr. Moreover, 
if codimension of is bigger than n, we can enlarge ^5^0 to a n-co dimensional foliation 
^, still horizontal with respect to the map tt giving rise to a common time-dependent 
superposition rule for the family ([1]). □ 
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4 Lie families, quasi-Lie and Lie systems 



This section is devoted to recalling the theories of quasi-Lie schemes and Lie systems 
needed to investigate the relations among these theories and Lie families. A full detailed 
report on these topics can be found in [71 [T7]. 

The theory of quasi-Lie schemes provides various results on the transformation prop- 
erties of time-dependent vector fields by a certain kind of time-dependent changes of 
variables associated with generalized flows. 

Each time-dependent vector field X gives rise to a generalized flow , i.e., a map 
g^ : (t, x) G M X — )■ gf{x) = g^{t, x) G M"' (more precisely, defined in a neighborhood 
of {0} X in M X M"), with g^ = Hku, such that the curve -f^^{t) = (xq) is the 
integral curve of the time-dependent vector field X starting from the point xq G M", i.e., 

7^0 = ^(^'7foW) and (0) = 9oM = Xo- 

Denote with Xt(M"') the set of all time-dependent vector fields on M". Each generalized 
flow h acts on the set of time-dependent vector fields Xj(M"') transforming each time- 
dependent vector field X G into a new one, h^X, with the generalized flow of the 
form g^*^ = h o g^ . In terms of autonomizations we can write ([T71 Theorem 3]) 

h^X = hjfX , 

where h is the natural autonomization of the generalized flow /i to a (local) diffeomorphism 
of M X M", h{t,x) = (t, ht{x)) and where /i* is the standard action of the diffeomorphism 
h on vector fields. This, in turn, implies that 

[7^,7^] = /i,[X,F]. (9) 

Let \^ be a finite-dimensional vector space of vector fields on M". We denote with 
V{M.) the set of time- dependent vector fields X G j£((]R") with values in V, that is, those 
time-dependent vector fields X such that, for every t G M, the vector field Xt{x) belongs to 
V. In terms of the introduced terminology and notation. Lie Theorem, whose statement 
can be found for instance in [U [7], can be reformulated as follows. 

Proposition 11. (Lie Theorem) A non- autonomous system X is a Lie system on M" 
if and only if there exists a finite- dimensional Lie algebra of vector fields Vq C j£(]R") such 
thatX G Vo{R). 

Definition 12. A quasi-Lie scheme S{W, V) on the manifold M consists of two finite- 
dimensional vector spaces of vector fields W,V C X{M) such that 

• W is a linear subspace of V. 

• W is a Lie algebra of vector fields, that is, [W, W] C W. 

• W normalizes V, i.e., [W,V] C V. 

It has been proved in [l7j that given a quasi-Lie scheme S{W,V), the space V{R) 
is stable under the action of the infinite-dimensional group Q{W) of generalized flows of 
vector fields in W{R), i.e., g^X G V{R), for every X G V(R) and g G g{W). 
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Definition 13. Given a quasi-Lie scheme S{W, V), we say that a time-dependent vector 
field X G V{M.) is a quasi-Lie system with respect to this scheme, if there exist a gener- 
ahzed flow g G G(W) and a Lie algebra of vector flelds Vo CV, such that g^X G Vo(M). 

As for each Lie system X there exists a Lie algebra of vector flelds Vq such that 
X G Vo(M), it is obvious that X G S{Vo, Vq) and, consequently, every Lie system is also a 
quasi-Lie system. 

From now on, given a quasi-Lie scheme S(W,V), a generalized flow g G Q(W), and 
a Lie algebra of vector flelds Vq C V, we denote with Sg(W, V; Vq) the set of quasi-Lie 
systems of the scheme S(W, V) such that g^X G Vo(M). 

Proposition 14. The family of quasi-Lie systems Sg{W, V; Vq) is a Lie family admitting 
the common time- dependent superposition function of the form 

^g{t,x^i), . . . ,x^rn),k) = g^^ o ^ [gt{x(^i)), . . . ,gt{x(^^)),k) , (10) 

for any time-independent superposition function $ associated with the Lie algebra of vector 
fields Vq by Lie Theorem. 

Proof. Let Zi, . . . , Zr he a basis in Vq. Since Vq is closed with respect to the Lie bracket, 

r 

[Zj,Zk\=^c,jkiZi (11) 

1=0 

for some constants Cjku j, k,l = 1, . . . ,r. For any Y G Sg{W, V; Vq), there exist functions 
bj such that 

r 

{g^Y),{x) = J2bAt)ZA^)- 

i=i 

Consequently, for the autonomization we can write 

r 

'^{t,x) = J2bAt)Zj{t,x), (12) 

j=0 

where we put Zq = (thus Zq = d/dt) and 6o(^) = 1 — Sj=i ^i(^)- 

Note that, as Zk are time- independent, the autonomizations Zk, k = 0, . . . ,r, form a 
Lie algebra: 

r 

[Zj, Zk] = Cjfc/Z/ , (13) 

1=0 

where cqm = cjqi = and cjko = — X^Li ^jkh for j. A; = 1, . . . , r. Hence, according to ([9]), 
the autonomizations = 5'^^(^fc) are also closed with respect to the bracket, 

r 

[Z'^,Z',] = Y,c,mZI, (14) 

1=0 
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and, in view of ( IT2l) . the autonomization of any Y G Sg{W, V; Vq) can be written in the 
form 

r 

j=0 

This means, in view of Theorem [T0| that Sg{W, V; Vq) is a Lie family. The form ( !T0|) can 
be now easily derived (see [T71 Theorem 4]). □ 

In view of the above proposition, every quasi-Lie system and, consequently, every Lie 
system can be included in a Lie family satisfying Theorem [101 This fact justifies once 
more calling this theorem Generalized Lie Theorem. 



5 Applications. 

In this section we will apply common time-dependent superposition rules for studying 
some first- and second-order differential equations. In this way, we will show how that 
common time-dependent superposition rules can be used to analyze equations which can- 
not be studied by means of the usual theory of Lie systems. Additionally, some new 
results for the study of Abel and Milne-Pinney equations are provided. 

5.1 A time-dependent superposition rule for Abel equations 

We illustrate here our theory by deriving a common time-dependent superposition rule 
for a Lie family of Abel equations whose elements do not admit a standard superposition 
rule except for a few particular instances. In this way, we single out that our theory 
provides new tools for investigating solutions of non-autonomous systems of differential 
equations than cannot be analyzed by means of the theory of Lie systems. 

With this aim, we analyze the so-called Abel equations of the first-type [201 121], i.e., 
the differential equations of the form 

doc 

— = ao{t) + ai{t)x + a2{t)x'^ + 03(^)2;^ (15) 

with 03 (t) 7^ 0. Abel equations appear in the analysis of several cosmological models 
[221 1231 111 and other different fields in Physics [2S1 il 123 I2S1 12S1 EO] . Additionally, the 
study of integrability conditions for Abel equations is a research topic of current interest in 
Mathematics and multiple studies have been carried out in order to analyze the properties 
of the solutions of these equations [211 EH [321 E31 El]- 

Note that, apart from its inherent mathematical interest, the knowledge of particular 
solutions of Abel equations allows us to study the properties of those physical systems that 
such equations describe. Thus, the expressions enabling us to obtain easily new solutions 
of Abel equations by means of several particular ones, like common time-dependent su- 
perposition rules, are interesting to study the solutions of these equations and, therefore, 
their related physical systems. 
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Unfortunately, all the expressions describing the general solution of Abel equations 
presently known can only be applied to study autonomous instances and, moreover, they 
depend on families of particular conditions satisfying certain extra conditions, see [311 
[32] . Taking this into account, common time- dependent superposition rules represent an 
improvement with respect to these previous expressions, as they permit one to treat 
non-autonomous Abel equations and they do not require the use of particular solutions 
satisfying additional conditions. 

Recall that, according to Theorem [TOl the existence of a common time-dependent 
superposition rule for a family of time-dependent vector fields ([2]) requires the existence 
of a system of generators, i.e., a certain set of time- dependent vector fields, Xi, . . . ^X^, 
satisfying relations (jT]). Conversely, given such a set, the family of time-dependent vector 
fields Y whose autonomizations can be written in the form 

r r 

= ^6,(t)X,(t,x), 5^6,(t) = 1, 

admits a common time-dependent superposition rule and becomes a Lie family. 

Consequently, a Lie family of Abel equations can be determined, for instance, by 
finding two time- dependent vector fields of the form 



d 

Xi{t,x) = {bo{t) + bi{t)x + b2{t)x^ + h{t)x') — , 

X2(t, x) = %it) + b[it)x + b'^{t)x^ + b'.it)^')-^, b'.it) ^ 0, 



(16) 



such that 

[Xi,X2] =2(X2-Xi). (17) 

Let us analyze the existence of such two time-dependent vector fields Xx and X2 
holding relation (ITTI) . In coordinates, the Lie bracket [Xi,X2] reads 

[(&3&2 - 62&3)x' + (2(6'3&1 - &3&'l) - &3 + K)^'' + (-3(6[,63 - 60^3) + " &26'i) 

- 62 + K)^^ + (-26062 + 2606^ - 61 + K)x - b'obi + 606; - 60 + b'o]-^. 

ox 

Hence, in order to satisfy condition ( ITTll . 6362 — 6363 = 0, e.g. we may fix 62 = 63 = 0. 
Additionally, for the sake of simplicity, we assume 63 = 1. In this case, the previous 
expression takes the form 

[26ix3 + (36o + 6'26i + b'2)x^ + (26o6'2 - 61 + b[)x - 6[,6i + 606^ - 60 + 6[,]^, 

and, taking into account the values chosen for 62, 63 and 63, assumption (ITTI) yields 61 = 1 
and 

6^ = 36o + 6^, 
2(6; - 1) = 2606^ + 6'i, 
2(&o - ^0) = -K + bob[ -bo + b',. 
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As the above system has more variables than equations, we can try to fix some values of 
the variables in order to simplify it and obtain a particular solution. In this way, taking 
bo(t) = t, the above system reads 

K = b2 - 3t, 
b[ = 2{b[ - 1) - 2tb'„ 
b'o = 2% -t) + b'Q- tb[ + 1. 

This system is integrable by quadratures and it can be verified that it admits the particular 
solution 

b'^{t) = 3{l + t), b[{t) = 3(1 + tf + 1, b'^it) = {1 + tf + t. 
Summing up, we have proved that the time-dependent vector fields 




((1 + tf + t+ (3(1 + + l)x + 3(1 + t)x^ + x^)^. 

ox 



(18) 



satisfy (fTTll and, therefore, the family of time-dependent vector fields Y'b(f)(t, x) = (1 — 
b{t))Xi{x) + b{t)X2{x) is a Lie family. The corresponding family of Abel equations is 

fix 

— = ^t + x) + b{t){l+t + xf. (19) 

According to the results proved in Section [3l in order to determine a common time- 
dependent superposition rule for the above Lie family we have to determine a first-integral 
for the vector fields of the distribution V spanned by the time-prolongations Xi and X2 
on M X ]R"('"+i) for a certain m so that the time-prolongations of Xi and X2 to M x M"*^ 
were linearly independent at a generic point. Taking into account expressions (fT8l) . the 
prolongations of the vector fields Xi and X2 to M x are linearly independent at a 
generic point and, in view of ( !T7|) . the time-prolongations Xi and X2 to M x span an 
involutive generalized distribution V with leaves of dimension two in a dense subset of 
M X M^. Finally, a first-integral for the vector fields in the distribution V will provide us 
a common time-dependent superposition rule for the Lie family fll9p . 

Since, in view of (fTTI) . the vector fields Xi and X2 span the distribution V, a function 
G : M X — )■ M is a first-integral of the vector fields of the distribution V if and only if 
G is a first-integral of Xi and Xi — X2, i.e. XiG = {X2 — Xi)G = 0. 

The condition XiG = reads 

dG , ,dG , ,dG ^ 

ar + ^^ + ^°^^ + ^^ + ^^)^ = °' 



and, using the method of characteristics ^35j, we note that the curves on which G is 
constant, the so-called characteristics, are solutions of the system 

, dxQ dxi dxi . ^ ^ 

dt = = =^ —— = t + Xi, 2 = 0, 1, 

t + xq t + Xi dt 
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i.e., Xi{t) = ^iC^ — t ~ 1, with i = 0,1. These solutions are determined by the imphcit 
equations ,^0 = e~*(a;o + t + l) and ^1 = e~^{xi + t + 1), with ^o, (,1 G M. Therefore, there 
exists a function G2 : — )■ M such that G{t, Xo,Xi) = G2{C,o,C,i). In other words, each 
first-integral G of Xi depends only on ^0 and ^1. 

Taking into account the previous fact, we look for first-integrals of the vector field 
X2 — X1 being also first-integrals of Xi, that is, for solutions of the equation (X2 — Xi)^ = 
with G depending on ^0 and ^1. Using the expression of X2 — Xi in the system of 
coordinates {t,^o,^i}, "we get that 

^dG ^dG _ ^dG2 ^dG2 

and, applying again the method of characteristics, we obtain that there exists a function 
Ga : M ^ M such that G{t, xq, Xi) = ^2(^0, 6) = G'3(A), where A = e^\{xo + 1 + l)-^ - 
(xi + t + 1)^^). Finally, using this first-integral, we get that the common time-dependent 
superposition rule for the Lie family (IT^ reads 

k = e2*((xo + t + _ (a;^ + t + i)-2)^ 

with k being a real constant. Therefore, given any particular solution Xi{t) of a particular 
instance of the family of first-order Abel equations fl2Tl) . the general solution, of this 
instance is 

x{t) = ((xi(t)+t + l)-' + A;e"'*)"'^'-t-l. 

Note that our previous procedure can be straightforwardly generalized to derive com- 
mon time-dependent superposition rules for generalized Abel equations [36], i.e., the dif- 
ferential equations of the form 

doc 

— = ao{t) + ai(t)x + a2{t)x^ + . . . + a„(t)x", n > 3. 

(JjV 

Actually, their study can be approached by analyzing the existence of two vector fields of 
the form 

Fi(t, x) = iPoit) + hi{t)x + ... + 6„(t)x")^, 

ox 

Y2{t,x) = {b',{t) + b[{t)x + ... + Kit)xn^, Kit) ^ 0, 

satisfying the relation [Yi,y2] = 2(1^2 — ^1) and following a procedure similar to the one 
developed above. 



5.2 Lie families and second-order differential equations 

Common time-dependent superposition rules describe solutions of non-autonomous sys- 
tems of first-order differential equations. Nevertheless, we shall now illustrate how this 
new kind of superposition rules can be applied to analyze also families of second-order 
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differential equations. More specifically, we shall derive a common time- dependent su- 
perposition rule in order to express the general solution of any instance of a family of 
Milne-Pinney equations [311 EZl EH] in terms of each generic pair of particular solutions, 
two constants, and the time. In this way, we provide a generalization to the setting of 
dissipative Milne-Pinney equations of the expression previously derived to analyze the 
solutions of Milne-Pinney equations in [TT] . 

Consider the family of dissipative Milne-Pinney equations [STJ |38l |39l HQ] of the form 

X = -Fx + uj'^x + e"2■^x"^ (20) 

with a fixed time-dependent function F = F{t), and parametrized by an arbitrary time- 
dependent function u = u{t). The physical motivation for the study of dissipative Milne- 
Pinney equations comes from its appearance in dissipative quantum mechanics [HI 
HU m], where, for instance, their solutions are used to obtain Gaussian solutions of 
non-conservative time- dependent quantum oscillators [l3]. Moreover, the mathematical 
properties of the solutions of dissipative Milne-Pinney equations have been studied by 
several authors from different points of view as well as for different purposes [TT | IT7 t [T51I571 
[3HI SSI SSI SZ] • As relevant instances, consider the works [TSl EZj which outline the state- 
of-the-art of the investigation of dissipative and non-dissipative Milne-Pinney equations. 
One of the main achievements on this topic (see [371 Corollary 5]) is concerned with an 
expression describing the general solution of a particular class of these equations in terms 
of a pair of generic particular solutions of a second-order linear differential equations and 
two constants. Recently the theory of quasi-Lie schemes and the theory of Lie systems 
enabled us to recover this latter result and other new ones from a geometric point of view 

[iniin]. 

Note that introducing a new variable v = x, we transform the family (120|) of second- 
order differential equations into a family of first-order ones 

X = V, 

7^ , 2 , -2F -S ^'^^^ 

V = —r V + 00 X + e x , 

whose dynamics is described by the following family of time- dependent vector fields on 
TM parametrized by u, 

= i-Fv + e-2^x-=^ + u\) — + V—, ueA = C^{t). 
V / ov ox 

Let us show that the above family is a Lie family whose common superposition rule can 
be used to analyze the solutions of the family (!20l) . 

In view of Theorem [101 if the family of systems related to the above family of time- 
dependent vector fields is a Lie family, that is, it admits a common time-dependent 
superposition rule in terms of m particular solutions, then the family of vector fields on 
M X R"(™+i) given by 



K., [VL, K.'], [K;, [K.', K."]], [K;, [K.', [K.", X.'"]]] . . . , a;, a;', uj'\ uj"\ . . . G A, (22) 
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spans an involutive generalized distribution with leaves of rank r < n ■ m + 1. 

^ Note that the distribution spanned by all is generated by the vector fields Yi and 

Y2, with 



-Fv + e'^^x'^ + x] ^ + v^, ¥2 



d_ 

dv 



d_ 

dx ' 



-Fv + e 



ov ox 



since = {1 — u'^)Y2 + ijJ^Yi. It is easy to see that the prolongation [Yi, Y2] is not spanned 
by Yi and Y2 and, so that we have to include the prolongation Y3 = [Yi, Y2] to the picture, 
where 

d ■ d 

I3 = X— (f + xF)- 



dx 



dv' 



In the case m = 0, the distribution spanned by the vector fields, Yi-, Yj,, I3, does not admit 
a non-trivial first-integral. In the case m > 0, the vector fields Yi, Y2, Y3 do not span all 
the elements of family fl22|) and we need to add to them the prolongation Y4 = [Yi, [Yi, 1^2]], 
with 

F4 = (2w + xF)— + {2e-'^^x-^ -2x- F{v + xF) - xF) — . 

ox ov 



The vector fields Yi, Y2, Y3, Y4 satisfy the commutation relations 



Yi,Y2 



Y3, 
n, 

{4 + F^ + 2F)% - {FF + F){Yi - Y2), 
2{Y, + 

{2 + F^ + 2F)% - {FF + F){Yi - 
-2% - 2(ri - F2)(4 + + 2F). 



Consequently, the vector fields ^1,^25^3,^4 span the vector fields of the family 
Adding Yi to each prolongation of the previous set, that is, considering the vector fields 
X\ = Yi, X2 =JY2,^X3^= Yi + Y3, and X4 = Yi + Y4, we get that the family of time- 
prolongations, Xi, X2, X3, X4, which spans the vector fields of the family fl22l) . The com- 
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mutation relations among them read 



Xi, X4 
X2, X3 
X2, X4 



X3 — Xi , 
X4 — Xi , 

-(FF + F + 4 + ^ 2F)Xi + (FF + F)X2 + (4 + F^ ^ 2F)X3, 
2X1 — 2X2 — X3 + X4, 

-(1 + F^ ^ 2F + FF + F)Xi + (FF + F)X2 + (1 + F^ + 2F)X3, 
-3X4 + (4 + F^ ^ 2F)X3 + (8 + F + FF + + 4F)X2+ 
SF^ -6F-FF 



+ (-9 - 3F' -6F-FF- F)Xi. 

As a consequence of Lemma [H we get that the vector fields Xi , X2 X3 and X4 close on 
the same commutation relations as the vector fields Xi, X2, X3, X4. Hence, in view of 
Theorem [T0| the family (1211) is a Lie family and the knowledge of non-trivial first-integrals 
of the vector fields of the distribution V spanned by Xi, X2, X3, X4 provides us with a 
common time- dependent superposition rule. 

Note that, as the vector fields Xi, Xi — X2 and their Lie brackets span the whole 
distribution V, a function G : M x TM^ — i- M is a first-integral for the vector fields of the 
distribution V if and only if it is a first-integral for the vector fields Xi and X2 — Xi. 
Therefore, we can reduce the problem of finding first-integrals for the vector fields of the 
distribution V to finding common first-integrals G for the vector fields Xi and Xi — X2. 

Let us analyze the implications of G being a first-integral of the vector field 

2 

Xi-X2 = 5^x, 



i=0 



d_ 

' dvi 



The characteristics of the above vector field are the solutions of the system 
dvQ dvi dv2 

= = , dxQ = 0, dxi = 0, dx2 = 0, dt = 0, 

Xo Xi X2 

that is, the solutions are curves in MxTM^ of the form s (t, Xq, Xi, X2, fo(s), f ^2(5)), 
with ^02 = XoV2{s) —X2Vo{s) and ^12 = XiV2{s) —X2Vi{s) for two real constants ^02 and ^12. 
Thus, there exists a function (72 : — )• M such that G{p) = G2(t, Xq, Xi,X2, C,o2, ^12), with 
p G M X TR^, ^02 = XqV2 — X2V0, and ^12 = X1V2 — V1X2. In other words, G is a function 
of t, a;o, Xi, 0:2, ^02, ^12- _ 

The function G also satisfies the condition XiG = which, in terms of the coordinate 
system {t, xq, Xi, X2, ^02^12, V2}, reads 



X G=—+ ^^^^^ ~ 
dt X2 



^02) dG ^ {x^V2 - 62) dG ^ dG 

H ~ r V2- 



dxo 



X2 



dxi 



' dxo 



X2 Xi 



X, 



X 



2/ J 



dG 
d^2 



02 



-2F 



f X2 




\xl 


xl) _ 



dG 



0. 
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That is, defining the vector fields 



+ 



dt X2 dxi X2 dxQ 



+ 



-F^02 - e- 



2F I "^2 ^0 



Xr, 



x-^ 



-FCi2 - e- 
d 



f X2 


_Xi\ 


d 


\x\ 


a). 


dii2 



d^G2 



xn d xi d d 

+ + 



X2 dXo X2 dxi 8X2 ' 

the condition XiG = imphes that H1G2 



V2B.2G2 = and, as G2 does not depend on 



V2., the function G must be simultaneously a first-integral for Si and S2, i.e., SiG = 
and 'E.2G — 0. 

Applying again the method of characteristics to the vector field S2, wc get that F can 
depend just on the variables t, ^02, ^2, A02 = xq/x2 and A12 = Xi/x2, that is, there exists 
a function G3 : — )■ M such that xq, a;i, 0:2, Vo, 'J^i, 1^2) = G2(t, a;o, 0:2, C02, G2) = 
Gi{t^ ^02, 62, A02, A12). 

We are left to check out the imphcations of the equation SiG = 0. Using the coordinate 
system {t, ^02, ^12, A02, A12, i'2, X2} and taking into account that G{t, xq, xi, X2, vq, vi, V2) = 
Gz{t, ^02- 'C125 A02, A12), the previous equation can be cast into the form SiG = ^TiGs + 
'^2Gz = 0, where 



1=0 



d 



dAi2 



- e-^ (A,-3 - A.2) 



d 



i2 



T2 = -Fe 



d 



12 



9e 



02 
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d d_ 
8^02 dt ' 



As G3 depends on the variables t, A02, A12, ^12, C02 only, we have that TiG — and 
T2G = 0. Repeating mutatis mutandis the previous procedures in order to determine the 
implications of being a first-integral of Ti and T2, we finally get that the first-integrals 
of the distribution V are functions of /i, I2 and /, with 



li = e'^^{xoVi - XiVo)'^ + 



and 



Xi) \Xq 



Xi 
X2 



i = l,2. 



I = e^^{xiV2 — X2V1Y + 
Defining V2 — e^V2, Vi — e^Vi and Vq — c^Vq, the above first-integrals read 



(XoVi - XiVo) + 



Xo 



1,2, 



18 



and 



[X1V2 



X2V1Y + 



X2J \Xi 



Note that these first-integrals have the same form as the ones considered in [TU] for k = 1. 
Therefore, we can apply the procedure done there to obtain that 



^0 = y hxj + k2xl + 2\J \i2[-{x\ + xf) + Ixjxl] 
with A12 being a function of the form 

kik2l + (-1 + kl + kl) 



(23) 



Xi2{ki,k2,I) 



/2_4 



and where the constants ki and k2 satisfy special conditions in order to ensure that xq is 
real [H]. 

Expression fl23|) permits us to determine the general solution, x{t), of any instance of 
family ( I20p in the form 

x{t) = \l kixlit) + k2xl{t) + 2^ \i2[-{x\{t) + x\{t)) + Ixl{t)xl{t) ] , (24) 



with 



Z^^^'\xi{t)x2{t) ~ X2{t)xi{t)f 



Xl{t) 
X2{t) 



X2{t) 
Xl{t) 



in terms of two of its particular solutions, xi(t), X2(t), its derivatives, the constants ki 
and ^2, and the time (included in the constant of the motion /). 

Note that the role of the constant / in expression ( I24l) differs from the roles carried 
out by kl and ^2- Indeed, the value of I is fixed by the particular solutions xi{t), X2{t) 
and its derivatives, while, for every pair of generic solutions xi{t) and 0:2 (t), the values of 
kl and ^2 range within certain intervals ensuring that x{t) is real. 

It is clear that the method illustrated here can also be applied to analyze solutions of 
any other family of second-order differential equations related to a Lie family by introduc- 
ing the new variable v = x. Additionally, it is worth noting that in the case F{t) = the 
family of dissipative Milne-Pinney equations (120|) reduces to a family of Milne-Pinney 
equations appearing broadly in the literature (see [H] and references therein), and the 
expression (12^ takes the form of the expression obtained in [11] for these equations. 



6 Conclusions and Outlook. 

We have proposed a generalization of Lie Theorem in order to characterize those 
families of non-autonomous systems of first-order ordinary differential equations, the so- 
called Lie families, that admit a common time-dependent superposition rule. We have 
studied the relations of quasi-Lie systems and Lie families. 
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In order to illustrate the usefulness of our achievements, we have derived common 
time-dependent superposition rules for studying dissipative Milne-Pinney equations and 
Abel equations. In the case of Abel equations, our result expresses the general solution of 
any particular instance of a Lie family of non-autonomous Abel equations in terms of each 
generic particular solution, a constant, and the time. In this way, we have initiated a new 
approach to study the solutions of these equations. Additionally, it is worth noting that 
the analyzed Lie family of Abel equations contains an autonomous instance admitting a 
special kind of superposition rule derived by Chiellini [31] . Unhke such a special superpo- 
sition rule, our common superposition rule does not require the use of particular solutions 
obeying any kind of extra condition and, therefore, it clearly represents an improvement 
with respect to Chiellini's technique. 

We have shown how common time- dependent superposition rules can be used to ana- 
lyze second-order differential equations by means of the study of a family of Milne-Pinney 
equations. More specifically, we have derived a common-superposition rule allowing us 
to obtain the general solution of any instance of such a family in terms of a generic pair 
of its particular solutions, their derivatives in terms of the time, and the time. Such 
an expression represents an interesting improvement with respect to previous results and 
methods, as it generalizes the superposition rule given in [11] for the usual Milne-Pinney 
equations to the dissipative case. 

We hope to get in the future new results on the theory of common time-dependent su- 
perposition rules and, additionally, to describe new applications, where our achievements 
can be used. 
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